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The application of Random Matrix Theory to the Dirac operator of QCD yields predictions for the probability 
distributions of the lowest eigenvalues. We measured Dirac operator spectra using massless overlap fermions in 
quenched QCD at topological charge u — 0, ±1 and ±2, and found agreement with those predictions — at least 
for the first non-zero eigenvalue — if the volume exceeds about (1.2 fm)^. 



In QCD with Nf massless quark flavors, chi- 
ral symmetry is assumed to be broken sponta- 
neously as SU{Nf)R^ SUiNf)L SU{Nf)H+L, 
which generates — 1 Goldstone bosons in the 
coset space SU{Nf ). They pick up a mass — 
so they can be identified with light mesons — in 
the case of small quark masses. Their dynamics 
is described by chiral perturbation theory as a 
low energy effective model of QCD. To the lowest 
order in the momenta and masses, the effective 
Lagrangian reads 



C[U] = ^Tr [d^Ud^U] 



;Tt[M{U + U^)] (1) 



where U{x) S SU{Nf ) and M is the (diagonal) 
quark mass matrix. The low energy constants iv 
and S appear here as free parameters. In a box 
of size V = the relations Ag^^ ^ L <C m^^ 
characterize the e-regime [1]. 

In that regime we consider the eigenvalues 
(EVs) iXn of the Dirac operator, resp. the di- 
mensionless variable z„ = A„yS. The spectral 
density is defined by p(A) = {J2n ^("^ ^ Ki))^ and 
in particular the microscopic spectral density pg 
is given by 



(2) 



In the sectors of topological charge ±v it can be 
decomposed as 



(3) 
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where p^^^ (z) is the probability distribution of the 
fc-th EV (excluding the zero EV). The application 
of Random Matrix Theory (RMT) to QCD pro- 
vides explicit predictions for the functions p'^j^^ {z) 
in the e-regime [2]. 

To test these predictions on the lattice, we sim- 
ulated quenched QCD with the Wilson gauge ac- 
tion and the Neuberger overlap operator Z?ov at 
zero quark mass [3]. At /3 = 6 resp. 5.85 the mass 
of the Wilson kernel, — was fixed by /i = 1.4 
resp. 1.6. In contrast to the continuum, where the 
EVs are imaginary, the EVs of Dov lie on a circle 
in C with center and radius /i. For comparison to 
the continuum predictions we mapped this circle 
stereographically onto the imaginary axis. 

The functions p^^\z) depend significantly on 

\v\] for instance, the peak of p^i\z) moves to 
larger values of z as \v\ increases. This ef- 
fect could not be seen convincingly, however, 
in simulations with staggered fermions (at least 
for the couplings investigated so far) [4]. On 
the other hand, previous studies with Ginsparg- 
Wilson fermions were consistent with the RMT 
predictions in QED2 [5], QED4 [6] and in QCD 
on a 4** lattice [7] . QCD studies on larger lattices 
were presented in Refs. [3,8]. 

Usually the test of RMT predictions was done 
by comparing the functions p^^\z) to the his- 
tograms of lattice data. To avoid the arbitrariness 
of the bin size in that procedure, we compared in- 
stead the cumulative density 
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In Fig. 1 we show the RMT curves for pi'^{z), 
z^l = 0, 1 and 2, compared to data from the fol- 
lowing lattices: 10^ at (3 — 5.85, 12* at /3 = 6 and 
8"* at /3 = 5.85. The corresponding physical vol- 
umes are « (1.23 fm)*, (1.12 fm)* and (0.98 fm)*. 
We see that the k = 1 EV tends to agree quite well 
with RMT, if the volume exceeds about (1.2 fm)*. 
It disagrees, however, if the physical volume is too 
small; then the data do not distinguish any more 
the different topological sectors, which is consis- 
tent with Ref. [7]. 
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Figure 1. The cumulative density of the first non- 
zero EV in the sectors = 0, 1 and 2. We 
compare the RMT predictions ( curves ) to the lat- 
tice data on volumes (1.23 fm)* (iop), (1.12 fm)* 
(center) and (0.98 fm)* (bottom). The quality of 
agreement increases with the volume. 



This comparison involves S as the only free 
parameter, which was optimized for agreement 
with RMT. In (1.23 fm)* and (1.12 fm)* we ob- 
tained very similar values, S = (253 MeV)'^ resp. 
(256 MeV)^. These values are consistent with 
the literature, although S is expected to diverge 
logarithmically at large V in quenched QCD [9]. 
Further volumes and higher statistics would be re- 
quired to verify this behavior; our statistics was 
too small to associate error bars with the above 
values for E. 

The range of z where the predictions are com- 
patible with the RMT curves raises gradually as 
the volume is enlarged. For V < {1 fm)'' this 
range is so short that it doesn't even capture the 
peak of the first non-zero EV, hence the RMT 
prediction is not applicable in that regime. How- 
ever, our largest volume is still not sufficient to 
capture well the k = 2 EV; we show results for 
our larger two volumes in Fig. 2. 
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Figure 2. The cumulative density of the second 
non-zero EV at = 0, 1, 2. We compare the 
RMT predictions (curves) to the lattice data in 
the volumes (1.23 fm)* (stars) and (1.12 fm)"* 
(squares). The agreement is worse than in Fig. 
1, since the relevant values of z are larger. 

Next we present results for the full spectral den- 
sity pi'^\z) of eq. (2) (summed over k). Fig. 3 
shows a histogram at = 1 on the 12"* lattice, 
which roughly follows the RMT prediction up to 
the second peak, before turning into the bulk be- 
havior p'f\z) — Cq -\- ClZ^ . 

We also studied the "unfolded level spacing dis- 
tribution" which was considered earlier on a 4* 
lattice [10]. There we found good agreement with 
the behavior specific for SU (3) in the fundamcn- 
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Figure 3. The spectral density of the low lying EV 
on a lattice of volume (1.12 fm)*. The histogram 
of the lattice data roughly follows the RMT pre- 
diction (dashed line) up to the second peak, before 
turning into the cubic bulk behavior (dotted line). 

tal representation that we were using. In this 
case, the resuhs agree well with the prediction on 
all the three lattice sizes that we considered; the 
lower limit for the volume does not seem to affect 
this quantity, see Fig. 4. However, this type of 
"unfolding" only keeps track of the relative order 
of the EV in the ensemble, hence it is not sensi- 
tive to the energy scale and not directly physical. 
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Figure 4. The unfolded level spacing distribution 
on a 8^ lattice at j3 — 5.85. The histogram follows 
closely the curve predicted for the gauge group 
SU{3) in the fundamental representation. 

In view of the future numerical exploration 
of the e-regime by means of Ginsparg- Wilson 
fcrmions, our main conclusions are 

• The volume should be larger than about 
(1.2 fm)^, which might just reflect that chi- 
ral perturbation theory is only applicable 
in the confinement phase. This observation 



can also be viewed as an empirical determi- 
nation of the so-called Thouless energy. 

• Once we are in the right regime, the RMT 
predictions for the low lying EVs are con- 
firmed in some range z — . . . Zcrit , where 
^ciit grows with the volume. 

• The latter means that particularly in the 
topologically neutral sector there is a non- 
negligible density of very small non-zero 
EVs. Such EVs are dangerous for the mea- 
surement of physical observables; they lead 
to the requirement of a huge statistics [11]. 

In the e-regime, observables tend to be strongly 
dependent [12], hence they should be mea- 
sured at fixed The last observation now im- 
plies that such measurements should beware of 

REFERENCES 

1. J. Gasser and H. Leutwyler, Phys. Lett. B188 
(1987) 477. 

2. P.H. Damgaard and S.M. Nishigaki, Nucl. 
Phys. B518 (1998) 495; Phys. Rev. D63 
(2001) 045012. T. Wilke, T. Guhr and T. 
Wettig, Phys. Rev. D57 (1998) 6486. 

3. W. Bietenholz, K. Jansen and S. Shcheredin, 
JHEP 07 (2003) 033. 

4. P.H. Damgaard, U.M. Heller, R. Niclasen and 
K. Rummukaincn, Phys. Rev. D61 (2000) 
014501. B.A. Berg, H. Markum, R. Pulhrsch 
and T. Wettig, Phys. Rev. D63 (2001) 014504. 

5. F. Farchioni, I. Hip, C.B. Lang and M. 
Wohlgenannt, Nucl. Phys. B549 (1999) 364. 

6. B.A. Berg et al., Phys. Lett. B514 (2001) 97. 

7. P. Hasenfratz et al., Nucl. Phys. B643 (2002) 
280. 

8. Contributions by T. Streuer and by P. Weisz 
to these proceedings. 

9. P.H. Damgaard, Phys. Lett. B608 (2001) 162. 

10. R.G. Edwards, U.M. Heller, J.E. Kiskis and 
R. Narayanan, Phys. Rev. Lett. 82 (1999) 
4188. 

11. P. Hernandez, K. Jansen and L. Lellouch, 
Phys. Lett. B469 (1999) 198. Contributions 
by T. Chiarappa and by K.-I. Nagai to these 
proceedings. 

12. H. Leutwyler and A. Smilga, Phys. Rev. D46 
(1992) 5607. 



